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CHARTS AND APPROXIMATE FORMULAS FOR THE ESTIMATION 
OF AEROELASTIC EFFECTS ON THE LATERAL CONTROL 
OF SWEPT AND UNSWEPT WINGS 
By Kenneth A. Foss and Franklin W. Diederich 


SUMMARY 


Charts and approximate formulas are presented for the estimation 
of static aeroelastic effects on the spanwise lift distribution, 
rolling-moment coefficient, and rate of roll due to the deflection of 
ailerons on swept and unswept wings at subsonic and supersonic speeds. 
Two types of stiffness distributions are considered, one which consists 
of a variation of stiffness with the fourth power of the chord and one 
which is based on an idealized constant-stress structure. Some design 
considerations brought out by the results of this paper are discussed. 
This paper treats the lateral-control case in a manner similar to that 
employed in NACA TN 2608 for the symmetric-flight case and is intended 
to be used in conjunction with NACA TN 2608 and the charts and formulas 
presented therein. 


INTRODUCTION 


The lateral control and maneuverability of a wing are important 
design considerations. These characteristics may be affected to a 
significant extent by aeroelastic action, particularly at high dynamic 
pressures and in the case of thin wings, swept wings, and wings designed 
for low wing loadings, because the operation of ailerons and spoilers 
usually creates aerodynamic forces which deform the wing. 

As a result of these deformations, the angles of attack along the 
span often change in such a manner as to produce lifts which oppose the 
rolling moment of the aileron or spoiler; furthermore, these lifts cause 
additional deformations which may again reduce the rolling moment, and 
so on, until equilibrium is reached. Wing flexibility may thus cause 
a serious loss in the control power; in fact, if the dynamic pressure 
of the air stream is sufficiently high, the aileron rolling moment may 
be completely nullified. The speed and dynamic pressure at this condi- 
tion are often referred to as the aileron reversal speed and reversal 
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dynamic pressure, because at higher dynamic pressures the controls would 
have to be reversed in order to roll the airplane. When wing flexibility 
causes a loss in lateral control, there is also usually a loss in the 
rolling maneuverability, which may be expressed as the wing-tip helix 
angle due to rolling and is affected by changes in both the control 
power and the damping in roll. These aeroelastic effects on the lateral 
control and maneuverability have to be taken into account in the design 
of a wing. 

Several methods are available for calculating these effects (refer- 
ence 1 , for instance), but since these effects depend on the structural 
characteristics of the wing, which are not accurately known in advance 
of its design, the relatively large amount of time required for even the 
most efficient of these methods militates against their use in connec- 
tion with preliminary design calculations. A need exists, therefore, 
for means of estimating some of the more important aeroelastic effects 
on lateral control quickly with an accuracy that is sufficient for pre- 
liminary design purposes. 

The related problem of estimating static aeroelastic effects on 
the magnitude and spanwise distribution of the lift in symmetric flight 
has been treated by the charts and approximate formulas presented in 
reference 2. The present paper consists of an extension of the analysis 
of reference 2 to the lateral-control case. Inasmuch as the static 
aeroelastic equations are linear, the results presented in the two 
papers may be superimposed. Included in the present paper are approxi- 
mate formulas for the estimation of the static aeroelastic effects on 
the spanwise lift distribution, rolling moment, and rate of roll due to 
aileron deflections on swept and unswept wings at subsonic and supersonic 
speeds. Also presented are summary charts which indicate whether a 
given design is likely to be affected by losses in lateral control. By 
means of these charts and approximate formulas as well as those of 
reference 2, the conventional procedure of designing a wing on the basis 
of certain strength criteria, checking it for aeroelastic phenomena, 
and then reinforcing it, when necessary, to meet the stiffness require- 
ments imposed by these phenomena can often be simplified greatly inasmuch 
as the effect of some of these phenomena can be estimated in advance of 
design. 

In order to keep the length of the paper to a minimum and to avoid 
a repetition of much of the material presented in reference 2, the 
present paper has been written in such a manner as to facilitate its 
joint use with reference 2 rather than to make it entirely self-contained. 
The use of the charts and approximate formulas presented herein is 
described and the limitations of the charts and the light they shed on 
some design problems are discussed. A numerical example is included to 
illustrate the use of the approximate formulas of this paper. A brief 
description of the calculations (based on references 1 and 2) on which 
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the charts and approximate formulas are based is contained in the 
appendix to supplement the more detailed derivations in references 1 
and 2. 


A 



a 


b 



Cl tte 

Cz 

V 

C *5 

C P 6 


d 

E 

e 

e i 


SYMBOLS 

aspect ratio (b^/S) 

swept-span aspect ratio (a/cos^a) 

location of section aerodynamic center measured from 
leading edge, fraction of chord 

wing span, inches 

chord (measured perpendicular to elastic axis), inches 

aileron chord, inches 

section lift-curve slope per radian 

effective wing lift-curve slope per radian 


rolling-moment coefficient 
damping-in-roll derivative 


rolling-moment coefficient due to aileron deflection 

location of chordwise center of pressure of lift 

produced by aileron deflection measured from leading 
edge, fraction of chord 


dimensionless sweep parameter 



tan^A 


•) 


Young*s modulus of elasticity, pounds per square inch 

location of elastic axis measured from leading edge, 
fraction of chord 


dimensionless moment arm of section lift about elastic 
axis (e - a) 



dimensionless moment arm of lift due to aileron deflec- 
tion about elastic axis (cp^ - e) 

root-stiffness function given in reference 2, 
equation (B25 ) 

allowable bending stress, pounds per square inch 

dimensionless function of distance along span used in 
approximate formulas for angle of attack due to 
aeroelastic action of ailerons 

modulus of rigidity, pounds per square inch 

wing thickness, inches 

section bending moment of inertia, inches 1 * 

mass moment of inertia of entire airplane about its 
longitudinal axis, inches 1 * 

mass moment of inertia of both wings about longitudinal 
axis of airplane, inches 1 * 

section twisting moment of inertia, inches 1 * 

dimensionless parameters used in approximate formulas 
for dimensionless dynamic pressures at aileron 
reversal given in table I 


dimensionless sweep parameter 


' 8 t (GJ) ] 
V e l c r (El), 


tan A ) ; the 


dimensionless parameter k/e is identical to k 
except that e^ is replaced by eg 


lift per unit distance along span, pounds per inch 

rolling moment on both wings, inch-pounds 

bending moment about an axis perpendicular to elastic 
axis, inch-pounds 


free-stream Mach number 


design load factor 
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P 

pb/2V 


rolling acceleration, radians per second per second 
wing-tip helix angle due to roll, radians 
dynamic pressure, pounds per square foot 

/ c L a e ;i.c r 2 St 2 cos A\ 

dimensionless dynamic pressure 1—2 2 )• 

\lkk (GJ) r / 

the dimensionless dynamic pressure eq* is identical 
to q* except that e]_ is replaced by 


C L CpS^^sin A 


dimensionless dynamic pressure 


a, 


(ei ) 


S 

s 


s* 

T 

t i 

V 

W 

w s 

y 

a 


“8 


total wing area, square inches 

distance along elastic axis measured from wing root, 
inches 

dimensionless distance along elastic axis (s/s-fc) 

accumulated torque about elastic axis, inch-pounds 

distributed torque due to inertia loading, inch-pounds 
per inch 

airspeed, feet per second 

design gross weight of airplane, pounds 

weight of primary structure of both wings, pounds 

lateral coordinate, inches 

angle of attack in planes parallel to plane of symmetry, 
radians 

angle of attack equivalent to unit aileron deflection, 
radians 


P angle of local dihedral, radians, or spanwise slope of 

normal displacement of elastic axis 

aileron deflection measured in planes parallel to air 
stream, radians 


5 
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A 

X 

P 

9 

0 ) 

A 

Subscripts.: 


moment-arm. ratio (e2/ei) 

structural-effectiveness factors defined in equa- 
tions (15) and table 1 of reference 2 

angle of sweepback at elastic axis 

wing taper ratio (c-t/c r ) 

density, slugs per cubic foot 

angle of structural twist in planes perpendicular to 
elastic axis, radians 


tip-stiffness -ratio parameter 


(El), 


(El), 


constant stress 


unit step function of distance along span 


D at divergence 

i at inboard end of aileron; inertia, in equation (36b) 

R at aileron reversal 

r at wing root 

s structural (due to structural deformation) 

t at wing tip • 

0 rigid wing (for q* — eq* = q = 0) 

USE OF THE CHARTS AND APPROXIMATE FORMULAS 


Summary of Method and Scope of the Calculations on Which the 
Charts and Approximate Formulas Are Based 


A brief description of the method and scope of the calculations is 
given here to indicate the limitations of the charts and approximate 
formulas. A detailed description of the method is given in the appendix. 
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Most of these calculations were performed by an extension, based 
on reference 1, of the matrix method of reference 2. This method con- 
sists in solving the differential equations descriptive of an elasti- 
cally deformed wing under aerodynamic loadings by numerical methods 
employing matrix techniques. Treated by this method were wings with 
three taper ratios (l, 0.5^ 0.2), one aileron span, two types of stiff- 
ness distributions, several values of the sweep parameters k and d 
which include sweptf orw&rd , unswept, and sweptback wings, and several 
values of the section moment -arm ratio e and the dynamic -pressure 
ratio q/q^. Calculated for each case were the dynamic pressure at 
aileron reversal and the changes in the spanwise lift distribution and 
rolling moment due to aileron deflection. For the constant -chord, 
constant— stiffness wings calculations were also performed by an exten- 
sion of the analytic method of reference 2, which consists in solving 
the differential equations exactly for these relatively simple cases. 
These calculations were made for two aileron spans and several values 
of the parameters k, d, e, and q/q E . 

Two important approximations that have been made in all calculations 
are as follows: 

(1) Aerodynamic induction effects at subsonic speeds are taken into 
account by an over-all reduction and, in most of the calculations, by 
rounding off the strip-theory loading at the tip and at the inboard end 
of the aileron in a manner described in the appendix; at supersonic 
speeds strip theory is used, with a small reduction at the tip in the 
matrix calculations 

(2) The rigid-body rotations imparted to a swept wing by its tri- 
angular root portion are taken into account by a suitable choice of an 
effective root 

These assumptions were made to reduce the number of independent 
parameters and to make the results more generally applicable. The most 
severe limitation on the use of the charts is probably imposed by the 
fact that calculations have been made for only two types of stiffnesses: 

(1) Stiffness distributions which vary as the fourth power of the 
chord, such as those of solid wings 

(2) Stiffness distributions associated with structures designed for 
a constant level of combined bending and torsion stress at every point 
on the span 

Except for solid wings and those with geometrically similar cross sec- 
tions, for which the stiffnesses vary as the fourth power of the chord, 
the stiffness distributions of any given wing depend on the detailed 
design of the wing and cannot be generalized easily. Consequently, 
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the constant-stress concept, which is discussed in detail in appendix B 
of reference 2, has been used to estimate stiffness distributions for 
some of the calculations. This concept constitutes an effort to relate 
the stiffness of a wing to its strength on the basis of the following 
assumptions: 

(1) The combined bending and torsional stresses are constant along 
the span 

(2) The bending and torsional stresses are combined in such a 
manner that the sum of the ratio of the actual to the allowable bending 
stress and the ratio of the actual to the allowable torsion stress is 
equal to unity when the margin of safety is zero 

(3) The structure is of the thin-skin, stringer-reinforced shell 
type and its main features do not vary along the span; for instance, 
the number of spars and their chordwise locations are constant along 
the span 

( 4 ) At the design condition the spanwise distribution of the applied 
loading is proportional to the chord 


Selection of Parameters 

Geometric parameters.- The geometric parameters used in the analy- 
sis are defined in figure 1. The location of the effective root indi- 
cated in this figure is discussed in reference 2. In the present paper 
the angle of aileron deflection 6 is defined as being measured in 
planes parallel to the air stream. This angle is equal to the product 
of the angle of rotation of the aileron about the hinge axis and the 
cosine of the sweep angle of the hinge axis. 

Although most of the charts and approximate formulas are based on 
a half-span outboard aileron (s* i = 0.5), the results of the analysis 
of the uniform wings with full-span ailerons (s*^ = 0) and "tip ailerons" 
(s*}— ■ »i) indicate that, except for the angle-of-attack distributions, 
the results based on a half-span aileron may be expected to be valid 
for outboard ailerons having spans which differ considerably from 
one -half. 

Aerodynamic parameters . - The aerodynamic parameters which enter the 
analysis are the effective wing lift-curve slope, the location of the 
wing aerodynamic center, the location of the chordwise center of pressure 
due to aileron deflection, and the angle of attack equivalent to unit 
aileron deflection. An effective wing lift-curve slope Cj^, applicable 

to basic lift distributions due to built-in twist, aileron deflection, 
roll, or aeroelastic twist, is approximately given by the relation 
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C ^a 


A cos A 


A + 4 — 2i cos A 

2rt 


at subsonic speeds, where c ^ is approximately given by 


( 1 ) 


c I a - ~j= (2) 

yl - Mq^cos^A 

The basis of equation (l) is explained in reference 2. At supersonic 
speeds the effective wing lift-curve slope is approximately 


4 cos A 


Ci^ = 

^Iq^cos^A - l 


(3) 


provided M q is greater than l/cos A. If Mq is greater than 1 but 
less than and not too close to l/cos A, equations (l) and ( 2 ) may be 
used in the absence of better information; however, the results obtained 
for this range of Mach numbers should be used with caution. 

The lift-curve slopes given by equations (l) and ( 3 ) should not be 
confused with the rigid-wing lift-curve slope or the damping- in-roll 
derivative; they are merely effective values suitable for aeroelastic 
calculations. Values of the rigid-wing lift-curve slope and the 
damping-in-roll derivative can be used in conjunction with the methods 
of the present paper and of reference 2 to obtain the values of the 
flexible-wing lift-curve slope and the damping- in-roll derivative because 
means are presented herein and in reference 2 for estimating the ratios 
of the flexible-wing to the rigid-wing values. For this purpose any 
experimental information concerning the rigid-wing values can be used* 
if none is available references 3, 4, and 5 may be used at supersonic’ 
speeds and reference 6 , at subsonic speeds. 

The local aerodynamic centers are assumed to be at a constant frac- 
tion of the chord from the leading edge, so that they are all equal to 
the wing aerodynamic center as a fraction of the mean aerodynamic chord. 
The moment arm 62 is then given by the relation 


(4) 


e-L = e - a 
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The local centers of pressure of the lift due to aileron deflection 
are also assumed to he a constant fraction of the chord from the leading 
edge; and the moment arm &2 i- 3 then given by 

e 2 = cp 5 - e (5) 

Theoretical two-dimensional values of the parameter cpg are presented 
in figure 2 for both trailing-edge and leading-edge ailerons at subsonic 
and supersonic speeds. At subsonic speeds the effect of finite span is 
to shift the center of pressure rearward. An appropriate value for this 
rearward shift may be estimated from the following relation, based on 
lifting-line theory for unswept elliptic wings with full-span ailerons: 



where the subscripts II and III refer, respectively, to two- and three- 
dimensional values. The use of the swept-span aspect ratio in place of 
A should serve to extend this approximate relation to swept wings. 

Theoretical two-dimensional values of ag, the angle of attack, 
equivalent to unit aileron deflection, are also given in figure 2. At 
low aspect ratios the values of ag for subsonic speeds tend to be 
higher than these two-dimensional values; as the aspect ratio approaches 0, 
ag approaches 1, at least in the case of wings without reentrant trailing 
edges. Experimental values of both ag and cpg are preferable to . . 
theoretical values if they are available. For spoilers these values 
always have to be obtained experimentally. 

The effective lift -curve slope and the values of the section moment 
arms vary with the free-stream Mach number; hence, the appropriate values 
must be used at each flight condition for which aeroe las tic calculations . 
are made. 

The airspeeds at which the lateral-control aeroelastic phenomena 
are of interest enter the calculations in the form of the corresponding 
dynamic pressures. These dynamic pressures, in turn, are expressed in 
dimensionless form by means of the relations 


C L rt e l c r 2s t 2cos A 

q * = •_!_ 

q 144 


(GJ), 


( 6 ) 
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or 


€q* 


e2C r ^s^^cos A 
144 (GJ) r 


(7) 


q = 


q 

i44 


C T c_s + ^sin A 

S> t 


(El) r 


( 8 ) 


The ratios of these quantities. 


and 



s t (GJ) r 
(El) r 


tan A 


(9) 


k 

€ 


q _ s t (GJ)j 
eg* e 2 c r (El) 


tan A 


( 10 ) 


are independent of the dynamic pressure and are very viseful for ana- 
lyzing the aeroelastic behavior of swept wings. Two other dimension- 
less parameters, which are independent of the dynamic pressure, enter 
the calculations: 


e 



( 11 ) 


and 


d = 


(GJ) r 

(El) r 


tan^A 


( 12 ) 


Structural parameters . - For the purposes of an aeroelastic analy- 
sis the wing structure is characterized by the location of the elastic 
axis, the magnitude and distribution of the bending and torsional 
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stiffnesses El and GJ, and the magnitude of the rigid-body rotations 
imparted to the wing by its root (taken into account in this paper only 
by the location of an effective root). The selection of these struc- 
tural parameters is discussed fully in reference 2. 


Preliminary Survey of Loss in Lateral Control 

The information contained in some of the charts and approximate 
formulas presented in the following sections of this paper has been 
summarized in figures 3 to 5 for the purpose of ascertaining in advance 
of more detailed estimates, if desired, whether the aeroelastic phe- 
nomena considered in this paper are likely to affect the design of the 
wing structure. This preliminary survey is not essential to any of 
the further calculations but may show them to be unnecessary in some 
cases. These figures pertain to constant-stress wings with half-span 
outboard ailerons. 

The charts of figures 3(a), 4(a), and 5(a) pertain to wings of 
taper ratios 0.2, 0.5, and 1.0, respectively, with the moment arm 

equal to 0 (corresponding roughly to subsonic flow conditions and 
elastic-axis locations fairly far back on the wing). These figures 
show the dynamic-pressure parameter q* defined by either equation (6) 
or 


q* 


(1 + X) 

18432 


e^ cos A 


(13) 


plotted against the sweep parameter k defined by equation ( 9 ) or 


, 1 + X. G Ah . a / , 1 \ 

k Tjv. tan A (14) 

2 E ti 19 e x 


for several values of the aileron effectiveness parameter C ^ /C ^ 

Ql Bq 

and for two values of the stiffness parameter 

± - (EI) r / B l c r f 

k 2 (GJ) r \ 8 t / 


( 15 ) 
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or 


_d_ = 32 z ftl9 e l\ 2 11 6^ 7 2 

k 2 (1 + x ) 2 g\ a a / n 8 T i9 T i 15 T i b 2 


where F^, is a root-stiffness parameter defined and given in reference 2 
and where the structural-effectiveness factors r\ & , t) v t) 6 , t\j, 1 ) 3 , tj 9 , 

and n 19 are also defined in reference 2 . (if, at the time a pre- 
liminary survey of aeroelastic effects is to be made, no information 
whatever concerning the wing stiffness is available, equations ( 13 ), 

(14), and (l 6 ) may be used; if an estimate of the root stiffnesses (GJ) 
and (El) r is available, equations ( 6 ), ( 9 ), and (15) maybe used.) P 

The charts of figures 3(b), 4(b), and 5(b) are the same as those 
of figures 3(a), 4(a), and 5(a), respectively, except that these charts 
are for wings with the moment-arm ratio 6 equal to unity (correspondirg 
to subsonic conditions with the elastic axis fairly far forward on the 
wing). 


The charts of figures 3(c), 4(c), and 5(c) pertain to wings with 
the moment arm e 1 equal to zero (corresponding roughly to supersonic 
flow conditions). These figures show the dynamic-pressure parameter gcP 
defined by either equation (j) or 


6q* = 


(1 + X ) 2 COS A 

18432 G nW h r q 

■i£ T ^ F r q a 


(17) 


plotted against the sweep parameter k/e defined by equation ( 10 ) 


or 


k 1 + X G Aa 

^ = ~i^r 7 tanA 

'19 2 


( 18 ) 


for two values of the stiffness parameter 
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d (El)r/egc r \ ^ 

(k/e ) 2 ’ (GJ) r V~/ 

d _ 32 E / T1 l9 e 2 N \ 2 T *6 T) 7 2 

(k/e) 2 (l + k) 2 G\A A J W\ 3 \ 2 


(19) 


( 20 ) 


The various lines of the charts of figures 3 to 5 designate the 
conditions at which a wing designed on the basis of strength considera- 
tions alone is likely to encounter changes in aileron rolling moment by 
various amounts due to wing flexibility. Those lines for zero rolling 
moment also designate the conditions at aileron reversal. These charts 
should be used in conjunction with the preliminary survey charts in 
reference 2. After it has been found, through the use of the charts 
in reference 2, to what extent the wing design is affected by such aero- 
elastic phenomena as divergence and aerodynamic -center shift, the same 
procedure can be used with the preliminary survey charts in this paper 
to ascertain whether the wing design is likely to be affected by lateral- 
control difficulties. If these charts indicate the likelihood of signif- 
icant aeroelastic effects on the aileron rolling moment, further calcu- 
lations are desirable. The charts and approximate formulas of this 
paper may be used for the preliminary calculations; once the structure 
has been designed, more refined methods, such as that of reference 1, 
may be used. 


Calculation of the Aeroelastic Phenomena Related to Lateral Control 

Analysis of the many solutions for the aeroelastic phenomena con- 
sidered in this paper obtained by the methods given in the appendix 
shows that the data can be summarized by means of approximate formulas. 
These formulas involve the aerodynamic, geometric, and structural 
parameters of the wing through the dimensionless parameters k, k/e , 
e, and d (equations (9), (10), (ll), and ( 12 ) , respectively) and 
through a series of constants Kp to Ky. The constants are functions 
of the taper ratio, stiffness distribution, and aileron span and are 
given in table I. As in reference 2, the form of these approximate 
formulas has been guided by considerations based on an idealized semi- 
rigid wing, and the actual values of the constants K3 to Ky (Kp and 

K 2 have been given in reference 2) were obtained by fitting the solu- 
tion for the functions Bp to Bg defined in the appendix by equa- 
tions (A18) and (A27) to their approximate expressions, equations (A32). 
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Dynamic pressure at aileron reversal . - The solutions for the 
aileron reversal -speed obtained by the methods given in the appendix 
can be summarized by approximate formulas which give the dimensionless 
parameters q* R , (eq*) R , or q R - that is, the values of the parameters 

defined in equations (6), (7), and (8) which correspond to the value of 
the dynamic pressure at aileron reversal - in terms of the parameters k, 
k/e, e, and d defined by equations (9), (10), (11), and (12), 
respectively. ' ' 

♦ 

An approximate formula for q* R is 


K i( 1 - % j d ) 

1 + (K^ + KgK 3 d)« + Ktjd + % | d + Kyk 


( 21 ) 


For very small values of the section moment arm e^ and the 

resulting large values of the parameters e and k, the following 
alternative forms of equation (21) are more convenient to use: 


(eq*) R = 


Ki(i - k 3 I a) 


7 * % + W * *5 7 d + k 6 H d + K 7 7 


( 22 ) 


and 


q R " 


% (l - K 3 I d) 


; + («* + w)! + «5 f + *6 f I + k 7 


(23) 


When the angle of sweep is zero, equations (21) and (22) reduce, 
respectively, to 



K 


1 


1 + Kj^e 


( 24 ) 
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and 


(*q*) R = ^rr (25) 

K k + J 


and when the moment arm e-j_ is zero, as it may he in supersonic flow, 
equations (22) and (23) reduce to 


and 



M 1 ~ k 3 i d ) 

+ K 2 K 3 d + K ? J 


(26) 


- • - k 3 e d ) 

q R “ e 

h + f + *7 


(27) 


The constants to are given in table I for wings having half- 
span outboard ailerons and taper ratios of 0.2, 0.5, and 1.0, for the 
two different types of stiffness distributions.^ These constants were 
found from the results of the numerical matrix method derived in the 
appendix. Also given in table I are these constants for uniform wings 
having full-span ailerons and tip ailerons calculated from the results 
of the analytic integration method of the appendix. Since values of 
the constants K-j_ to K y are given for three aileron spans in the case 
of the uniform wing, they may be interpolated to yield values for other 
aileron spans. No calculations were made for other than half-span 
ailerons on tapered wings; nevertheless, as pointed out previously, 
calculated for half -span ailerons should be reasonably valid for out- 
board ailerons having spans which differ considerably from one -half. 

No aileron-reversal calculations have been made for swept wings with 
inboard ailerons. However, for unswept uniform wings the dynamic pres- 
sure at aileron reversal has been calculated for the limiting case of a 
wing with an inboard aileron of vanishingly small span by operational 
methods similar to those described in the appendix in connection with 



NACA TN 21 hi 


17 


the calculations for uniform wings. The value -of q* R obtained in 

this manner is shown as a function of e in figure 6. Also shown in 
figure 6 are the values of **R for full-span and for tip ailerons. 

For small values of e, such as are likely to be encountered in subsonic 
flight, there is little difference in the values of q*^ for the three 

aileron configurations, but at large values of 6, such as are likely 
to be encountered at high-supersonic speeds, there is some difference 
between them; the aileron reversal speed is highest for the tip aileron 
and lowest for the inboard aileron. However, these conclusions may 
not be valid for nonuniform or sweptback wings. 

With the values of q* R , (eq*) R , and q R given by equations (21), 

(22), and (23) and the definitions of these parameters given by equa- 
tions (6), (7), and (8), the values of q at aileron reversal may be 
determined. ■ If desired, the corresponding airspeed may be determined 
from the relation 

Vr \/| 

If both q R (as obtained from reference 2, for instance) 
are positive and q^ is greater than q R , there is no actual 

reversal speed; in fact, the aileron rolling moment will tend 
increase with dynamic pressure until divergence is reached. 

The value of q R calculated for any given value of q* R , 
or q^ depends on the value of the effective lift-curve slope Ci ^ 

and, hence, on the Mach number. As suggested in reference 1, the value 
of q R may be plotted against Mach number on log- log coordinates; if 
the straight lines of the actual dynamic pressure at several altitudes 
as functions of Mach number are drawn on the same plot, an intersection 
of the reversal line with one of the lines of actual dynamic pressure 
designates possible aileron reversal at that value of dynamic pressure, 
Mach number, and altitude. 

Spanwise angle-of-attack distribution . - In the appendix, an approxi- 
mate expression is determined for the change in angle of attack due to 
the deflection of ailerons on flexible wings. The ratio of the angle- 
of-attack distribution due to structural deformation a g to the 
effective angle of attack of an aileron agb is 


(28) 


and q R 
aileron 
to 

(eq*) R , 
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056 


q/q D (y + Kgt)^ - Kgk Af a 
_ q_ 1 - K 2 k 

q D 


(29) 


The functions f a and Af a , which depend on the spanwise coordinate s* 
are given in figure 7 for wings having half-span outboard ailerons, 
taper ratios of 0.2, 0.5, and 1.0, and the two different types of stiff- 
ness distribution. The value of . q^ required in equation (29) may be 
found from the approximate expression for q* D or q D given in refer- 
ence 2 as 


K, 


i* = 




(30) 


or 


Kl 



(31) 


When is very large, a more convenient form of equation (29) is 

as follows: 


eq . ( K 4 * *2 tK - K 2 I Af a 

a 5 5 1 - -i- K x 

q D 


(32) 


Spanwise lift distribution .- Within the limitations of the modified 
strip theory used in the analysis, the lift per inch of span is propor- 
tional to the local effective angle of attack, so that 



(33) 


where ccg/ccgS is obtained as indicated in the preceding section and 
J is a unit step function of the distance along the span defined by 


2a. ~ 0 ( when 3 < 3 i) 

i a = 1 (when s > 8^) 
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The expressions for angle -of -at tack deformations and lift distribu- 
tions (equations ( 29 ) and (33)) can be used only for half-span out- 
board ailerons or spoilers since the functions f a and Af a have been 

calculated only for this case. 

Rolling-moment coefficient. and rate of roll due to aileron 
deflection . - The rolling-moment coefficient due to a unit aileron 
deflection may be obtained in terms of its equivalent rigid-wing 

'value from the approximate formula 





(34) 


and the wing-tip helix angle due to an aileron deflection pb/2V, which 
is a measure of the rate of roll and the rolling maneuverability, may 
be obtained from 



(35) 


The manner in which the aileron-effectiveness parameter 


varies with dynamic pressure depends on the ratio of the reversal to 
the divergence dynamic pressure as may be seen from equation (3*0 and 
figure 8. When the aerodynamic center is ahead of the elastic axis, 
as is generally the case at subsonic speeds, i s positive and 

greater than one for swept forward wings, positive and less than one for 
unswept wings, and negative for sweptback wings. Thus, in general, 

increases with dynamic pressure until divergence is reached 


for swept forward wings; it decreases slowly at first and then more 
rapidly as reversal is approached for unswept wings; and it decreases 
rapidly at first and then more slowly for sweptback wings. The rapid 
decrease in rolling effectiveness for sweptback wings can be alleviated 
by the use of unconventional lateral-control devices which have their 
centers of pressure ahead of the elastic axis, such as leading-edge 
ailerons or spoilers. These devices may also serve to make the dynamic 
pressure at aileron reversal negative, so that there is no reversal of 
lateral control in the given speed range. The loss or gain of lateral 
control is then given by the part of figure 8 for negative values of 
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The analysis summarized by the approximate formulas ( 3*0 and ( 35 ) 
is based on rolling moments about the wing root instead of about the 
fuselage center line, partly to simplify the analysis and partly to 
avoid the introduction of the fuselage width as another independent 
parameter. However, these approximate formulas should be valid for 
obtaining rolling moments and rates of roll about the fuselage center 
line as well, because equations (3*0 and (35) are expressed as ratios 
of flexible-wing to rigid-wing values; that is, the ratio of the 
rolling moment about the fuselage center line to its rigid-wing value 
should be nearly the same as the ratio of the equivalent rolling 
moment about the wing root to its rigid-wing value when the ratio of 
the fuselage width to the wing span is small. 

, The rolling-moment coefficient and rates of roll given by equa- 
tions (3^) and (35) are functions of the aileron span inasmuch as q R 
is a function of the aileron span. The variation of C or pb/2V 

with aileron span can therefore be found only for uniform wings; how- 
ever, since the effect of aileron span on q R is not very great,, its 
effect on Cjg or pb/2V is not likely to be great. In the case of 

uniform wings, the rolling moment due to the deflection of an inboard 
aileron can be found by superposition because the aeroelastic equations 
upon which the results of this paper are based are linear; that is, the 
rolling-moment coefficient due to a 3 0-percent -span inboard aileron, 
for example, is equal to for a full-span aileron minus C jg for 

a 70-percent -span outboard aileron. 

Inertia effects .- In steady rolling flight no inertia effects are 
present which can affect the static aeroelastic problem except, possi- 
bly, for centrifugal forces on heavy underslung nacelles. The maximum 
value of pb/2V is therefore usually unaffected by inertia effects. 
However, in the equally important problem of initial rolling accelera- 
tion, which governs the time in which a given rolling velocity can be 
attained, inertia effects must usually be taken into account. In 
reference 2 the observation was made that in symmetric flight inertia 
effects are not as important as other static aeroelastic effects, 
except for flying wings, because the inertia forces are in about the 
same ratio to the aerodynamic forces as the wing weight is to the air- 
plane weight. By the same reasoning, inertia effects are almost always 
very important in getting into a roll because the inertia forces are 
then in about the same ratio to the aerodynamic forces as the moment 
of inertia of the wings about the longitudinal axis of the airplane is 
to the moment of inertia of the entire airplane about its longitudinal 
axis, a ratio which is usually close to 1. 

No charts are presented in this paper for these inertia effects 
because the manner in which mass is distributed varies so widely among , 



NACA TN 27^7 


21 


different airplanes that preparation of a generally applicable set of 
charts appears t,o be impractical at present. However, the procedure 
outlined in reference 2 for taking inertia effects into account in the 
calculation of quasi-static aeroelastic phenomena by means of the charts 
presented therein may be applied to the calculation of the inertia 
effects encountered in starting a roll. This procedure is described in 
the following paragraphs. 

For a given rolling acceleration p, the linear normal accelera- 
tion of an element of mass at a distance y from the center line of 
the airplane is py. From this linear acceleration and the known or 
estimated mass distribution of the wing, .the inertia load Z-^ per inch 
of span and the inertia torque t^ per inch of span can be calculated 
for any given normal, pitching, or rolling acceleration. Substitution 
of these loads and torques for the terms Z and Ze^c in equations (A3) 

or (A36) and equations (A2) or (A35) of reference 2, respectively, 
yields the values of the accumulated bending moments and torques in 
equations (A4) and (A5) or in equations (A37) and (A38) of reference 2. 
Equation (a 6) of reference 2, or the matrix equivalent of this equa- 
tion, then yields the angle -of-attack distribution due to the deforma- 
tions caused by the inertia effects associated with the given acceleration. 

This angle-of-attack distribution can be considered as a geometri- 
cal angle-of-attack distribution. For' the purpose of calculating the 
increment caused by aeroelastic action, this distribution can be approxi- 
mated by a linear-twist angle-of-attack distribution with a value at the 
wing tip which is such that the moment about the effective wing root of 
the area under the linear-twist distribution equals the moment of the 
area under the calculated angle-of-attack distribution due to inertia 
effects. (The moment, rather than the area, is suggested as a basis of 
correlation because the angles of attack near the wing tip are more 
important in aeroelastic phenomena than those at the wing root. ) The 
justification for this rather arbitrary approximation to the angle-of- 
attack distribution is that the correction to be applied as a result of 
aeroelastic action to the deformations due to inertia loads is usually 
small compared with these deformations. 

The angle of attack due to structural deformation a s associated ' 
with the linear-twist distribution can then be obtained from equa- 
tion (21) and figure 7 of reference 2 or, if X = 0, from figure 8 of 
reference 2. The lift distribution associated with the total angle-of- 
attack distribution due to the deformations caused by the inertia effects, 
including the increment in this angle-of-attack distribution produced by 
aeroelastic action, can then be found from equation (24b) of reference 2, 
in which' a,g, a^, an ^ Zq pertain to the calculated angle-of-attack 

distribution due to the inertia effects (not the linear approximation to 
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this distribution). This lift distribution can be integrated to obtain 
the rolling moment due to inertia effects, as modified by aeroelastic 
action. 


• The rolling moment calculated in this manner may then be combined 
with the rolling moment due to aileron or spoiler deflection, which 
may be calculated as indicated in the preceding section. If the con- 
tributions of the tail and the fuselage to the rolling moment are 
neglected. 


✓ 


- 1 *> - ^ c » 


5qSb + 
°s 



where p is the angular acceleration in roll, I is the mass moment 
of inertia of the entire airplane about its longitudinal axis, b w is 

the mass moment of inertia of both wings about the longitudinal axis 
of the airplane, and C j is the flexible-wing value of the rolling- 
's 

moment coefficient due to aileron deflection (which may be calculated 
in the manner described in the preceding section). The ratio (dL' w /dp) s 

is the rolling moment per unit roiling acceleration due to inertia 
effects, including aeroelastic effects, and is equal to -I w plus the 

rolling moment due to the lift which results from the deformations due 
to the inertia loads per unit angular acceleration in roll as well as 
from the aeroelastic deformations which accompany these inertia defor- 
mations; in other words, (dL’ w /<3p) is equal to -I w plus the rolling 

moment calculated as described in the preceding paragraphs for p = 1 . 
Then 


or 


C U &c ^ Sb 



( 36 a) 


5qSb 


P = 
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( 36 b) 
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where 



is a rolling-moment coefficient per unit aileron deflection which 
includes static aeroelastic effects, inertia effects, and the aero- 
elastic magnification of the inertia effects. This rolling^-moment 
coefficient is a truer index of the rate at which a roll can be 
initiated than Cjg . The initial rolling acceleration (disregarding 

unsteady-lift effects) can be calculated from equations ( 36 ). 


Illustrative Example 

The approximate formulas described in the preceding sections have 
been used to find the effects of aeroelastic ity on some lateral-control 
properties of the wing considered in the illustrative example of -refer- 
ence 2. The resulting calculations are an 'extension of those in refer- 
ence 2, and the additional parameters are presented in table II. 

The subsonic and supersonic values of the parameters k, e , and d 
were calculated from equations ( 9 ), (ll), and (12), respectively. With 
the appropriate values of the factors to K ^ interpolated from 

table I, the values of q*^ were calculated from equation (21) and are 
given in table II. From these values of q*-^, "the subsonic and super- 
sonic dynamic pressures at aileron reversal were found by means of 
equation (6). These values of q^ vary as the reciprocal of the 

effective lift-curve slope, if the corresponding values of e-^ and e 2 
are assumed to remain constant. 

In order to find the angle -of-attack distributions due to deflec- 
tions of the aileron from equation ( 29 ), the values of the functions f a 
and Af a were taken from figure 7(c).- The spanwise change in angle of 
attack is shown in the top plot of figure 9 for different values of the 
dynamic-pressure ratio. The rolling-moment coefficient and wing-tip / 
helix angle due to deflections of the aileron were calculated from 
equations (3*0 and (35) and were plotted in figure 9 as functions of 

q 

the dynamic-pressure ratio . 
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DISCUSSION 

Limitations of the Charts and Approximate Formulas 


The charts and the approximate formulas presented in this paper 
are subject to certain limitations as a result of the approximations 
made in the calculations on which they are based. These limitations 
are discussed fully in reference 2 and can be classified as restrictions 
on the plan form, on the speed regime, and on the wing structure. The 
limitations are given very briefly as follows: 

(1) The results obtainable by the use of the charts and approximate 
formulas are likely to be unsatisfactory for wings of very low aspect 
ratio, very large sweep, or zero taper ratio 

(2) The results are restricted to wings on which the spanwise lift 
distribution is roughly proportional to the chord and angle of attack 
and on which the section aerodynamic centers are at an approximately 
constant fraction of tne chord; these restrictions are most likely to 
be violated by wings flying at transonic speeds and by wings having 
concentrated sources of lift, such as nacelles and tip tanks 

(3) The results are somewhat restricted to wings with one of the 
two types of spanwise stiffness distributions used in the analysis, 
wings with no chordwise bending (relatively thick wings), and wings 
having an elastic axis at an approximately constant fraction of the 
chord. 

The manner in which the aeroelastic effects of aileron deflection 
are analyzed in the present paper imposes certain additional limitations 
that particularly affect the aileron geometry. In the analysis in the 
appendix, the spanwise lift distribution due to the deflection of an 
aileron was approximated by strip theory, an approximation which is 
probably less valid for aileron deflections than for geometric angles 
of attack. The assumption was also made that the centers of pressure 
due to aileron deflections are at a constant fraction of the wing chord; 
this assumption is also probably less valid than the assumption that the 
section aerodynamic centers of the wing are at a constant fraction of 
the chord. Since these assumptions are more nearly true for wings of 
high than those of low aspect ratio, these limitations serve, in effect, 
to restrict the applicability of the present paper to aspect ratios 
somewhat higher than those amenable -to the analyses of reference 2. 

The results of the present paper do not take into account explicitly 
any flexibility of the aileron itself because of the assumption that the 
angle between the aileron and wing is constant along the span of the 
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aileron. This assumption is almost universally made in analyzing the 
aeroelastic properties of ailerons and is justifiable because the net 
effect of the difference between the wing deformations and the aileron 
deformations on the over-all lift and moments appears to be negligible. 

As a result of the fact that the static aeroelastic phenomena 
associated with lateral control involve many more parameters than do 
those associated with symmetric flight, the coverage 'of the various 
parameters is not as complete as in reference 2. Specifically, angle- 
of-attack distributions have been presented only for uniform wings, all 
the charts and approximate formulas presented in this paper are 
restricted to outboard lateral-control devices (ailerons or spoilers) 
except for the uniform-wing case, and most of the calculations upon 
which these results are based were made for wings with half-span out- 
board ailerons (s* i = 0.5). However, the results of the analysis for 

the uniform wings with full-span ailerons (s^ = 0) and tip ailerons 
(s* i — >l) indicate that, except for the angle-of-attack distributions, 
the results based on a half-span aileron are approximately valid for 
outboard ailerons of spans differing considerably from one -half. 


Relation between Strength and Stiffness as 
Design Criteria 

The' relation between strength and stiffness as design criteria 
was discussed in reference 2. The preliminary survey charts in refer- 
ence 2 indicate the extent to which a wing designed on the basis of 
strength considerations alone is likely to be affected by aeroelastic 
phenomena and, consequently, indicate whether the wing has to be 
stiffened beyond the amount associated with the required strength. The 
preliminary survey charts of this paper- (figs. 3 to 5) serve the same 
purpose in regard to the aeroelastic effects on lateral control; 
furthermore, even though the charts of reference 2 may indicate that a 
particular wing is not significantly affected by the aeroelastic phe- 
nomena considered in that paper, this wing may still have to be 
stiffened because of an undesirably large loss in lateral control. 

As may be concluded from the survey charts of reference 2 and the 
present paper, as well as from the discussion contained in reference 2, 
the following wings designed on the basis of strength considerations 
are most likely to be subject to adverse aeroelastic effects on lateral 
control and rolling maneuverability: 

(1) Wings with a high design flying speed or dynamic pressure 

(2) Sweptback wings 
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(3) Thin wings 

(1<-) Wings designed for low wing loading 

(5) Unswept and moderately swept wings with an elastic axis 
relatively far forward on the chord or with the center of pressure of 
the lift produced by aileron deflections relatively far back on the 
chord as a result of the aileron configuration (small aileron chord or 
wing of low aspect ratio) or flight condition (supersonic speeds) 

(6) Wings with a relatively high lift-curve slope 


Structural Weight Associated with the 
Required Stiffness 

When a given wing has been shown to be subject to undesirably large 
aeroelastic effects by means of the charts of reference 2 and the present 
paper or by any other method, the problem arises how to distribute the 
additional required stiffness; that is, which spanwise distribution of 
structural material will alleviate the adverse aeroelastic phenomena to 
the desired extent with the minimum increase in structural weight. 

In order to shed some light on this problem, aeroelastic and weight 
calculations have been made for a family of somewhat arbitrarily selected 
stiffness distributions which differ from the distribution required . 
by the constant-stress criterion in a manner described in reference 2. 
These stiffness distributions are designated by the tip-stiffness -ratio 
parameter to, which is the ratio of the stiffness El or GJ at the 
wing tip to the corresponding stiffness of a constant-stress wing. The 
results of the lateral-control calculations for wings with taper 
ratio 0.5, with constant wing-thickness ratio h/c along the span, and 
with two of these stiffness distributions are included in table I and 
figure 7(b). The designation "excess strength" refers to the stiffness 
distribution increased over the constant-stress requirement to such an 
extent that the value of o> is 2.0. The results of the aeroelastic 
calculations for the stiffness distributions decreased below the constant- 
stress requirement to a value of co = 0.5 happen to be the same as the 
results for the constant-stress stiffness distributions for wings with 

(h/c) t 

varying wing-thickness ratio; that is, = 0.5. The structural 

(h/c) r 

weight considered in these calculations is that of the primary load- 
carrying structure; the remaining structure is assumed to be unchanged 
in the stiffening process. 
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The results of the weight calculations and aeroelastic calcula- 
tions in reference 2 indicated that the addition of stiffness in the 
outboard regions of a wing (large values of <x>) was more efficient, 
from weight considerations, in alleviating the aeroelastic effects 
considered in that paper than the addition of stiffness in the. inboard 
regions. This conclusion is corroborated, in essence, by the calcula- 
tions made for the aeroelastic phenomena considered in the present 
paper. Figure 10, which consists of a plot of the structural weights 
required for a given loss in lateral control at a constant value of 
dynamic pressure, indicates that the least weight is associated with 
values of the tip-stiffness-ratio parameter a) greater than 1, except 
for wings with values of the sweep parameters k or k/e equal to -8. 
These large negative values of k or k/e, however, pertain to wings 
that are (l) swept forward, (2) sweptback with the aerodynamic center 
behind the elastic axis (in the case of negative k), or ( 3 ) sweptback 
with the center of pressure due to aileron deflection ahead of the 
elastic axis (in the case of negative k/e, as it may be for spoilers 
or leading-edge ailerons). For swept forward wings, the aileron rolling 
moment usually increases rather than decreases with dynamic pressure; 
lateral control then does not impose any structural requirements so 
'that the curves do not have much significance for such wings. 


The Aeroisoclinic Wing 

It was shown in reference 2 that an over-all type of aeroisoclin- 
icism in which bending and torsion action tend to cancel for the wing 
as a whole can be achieved for the aeroelastic phenomena considered in 
that paper by a choice of a suitable ratio of the bending to the tor- 
sion stiffness or by a choice of the elastic-axis location, that is, by 
satisfying the relation 


s t (GJ)r 
e l C r (El) r 


tan A = 


K 2 


(37) 


(where Kg is given in table I). However, reference 2 indicates that 

even if the CQndit ions of equation (37) are achieved, there may still 
be great losses in lateral control and the wing may still be subject 
to adverse dynamic phenomena; in fact, the severity of adverse aero- 
elastic effects in the lateral control and of certain dynamic phenomena 
may be increased as a result of achieving aeroisoclinicism. 

The results of the present paper corroborate the conclusion con- 
cerning aeroelastic effects on lateral control. However, by suitable 
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additional modifications a wing which has been made aeroisoclinic can 
also be made to suffer no loss in lateral control due to aeroelastic 
action. As may be seen from figures 3 to 5 or from equation (3*0, the 
condition for no losses in lateral control is that be equal to q^; 

or, by setting equation (21) equal to 4 equation (30), that 


K^e + (Kg + K^)k + K^d + (K 3 



(38) 


of which the condition 


8 t (gJjr 
e 2 c r (El) r 


tan A = 


K2 


(39) 


usually is an approximate solution. Therefore, in order to satisfy 
both of the conditions specified in equations (37) and (38), the sec- 
tion moment arm e 2 must be nearly equal to -e-^ or, in other words, 
the center-of-pressure parameters a and cp^ must be nearly equal, 
a condition which can be satisfied by using full-chord ailerons (all- 
movable wing tips), or a combination of geared leading- and trailing- 
edge ailerons, for instance. However, as pointed out in reference 2, 
attempts at solving static aeroelastic problems by aiming at aeroiso- 
clinicism may tend to aggravate certain dynamic phenomena. The same 
statement must also be made concerning the foregoing methods of 
alleviating static aeroelastic effects on lateral control; these 
methods may, for instance, lead to flutter difficulties which may 
require excessive mass balancing of the control surfaces. 


Relation of Charts to Design Procedure 

The conventional procedure of designing a wing on the basis of 
strength requirements and later checking it for aeroelastic effects 
can be facilitated at several stages by using the methods described 
in the present paper and in reference 2. As pointed out in reference 2, 
for instance, the preliminary-survey charts presented therein can be 
used to establish some static aeroelastic characteristics that would 
be obtained in symmetric flight if the wing were designed for strength 
alone . . . 

If these characteristics are deemed satisfactory, the design can 
proceed on the basis of strength requirements alone. If, on the other 
hand, they are considered unsatisfactory, the wing must be stiffened. 
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The amount of additional material required can be est ima ted, as indicated 
in reference 2, by interpolating between the results presented therein 
for the constant -strength case, the "excess strength" case, and the 

(h/c ) t 

case TT = °*5. As previously mentioned, the additional structural 

0v<0 r 

material is usually most effective if distributed near the wing tip. 

Similarly, the preliminary-survey charts of the present paper can 
be used to ascertain whether the wing can meet lateral-control require- 
ments if designed for strength alone. If it must be stiffened to meet 
these requirements, the necessary amount of additional material can be 
estimated in the same manner as indicated in reference 2 for aeroelastic 
effects incurred in symmetric flight. 

Inasmuch as the charts of reference 2 and of the present paper 
pertain only to static aeroelastic phenomena, the problem remains of 
ascertaining in the preliminary design stage whether a wing designed 
for strength alone (or, for that matter, a wing designed both on the 
basis of strength requirements and of static aeroelastic considerations) 
is likely to experience flutter difficulties. However, flutter is a 
much more complicated phenomenon and depends on many more parameters 
than do static aeroelastic phenomena. Consequently, preparation of a 
.generally applicable set of charts appears impractical. Nonetheless, 
although phenomenologically or functionally flutter is not related to 
the static aeroelastic phenomena considered in reference 2 and the 
present paper, it is "mechanically" related by virtue of the fact that 
all these phenomena depend on the wing geometry and the wing stiffness 
(althqugh the aerodynamic parameters are different and flutter, unlike 
the static aeroelastic phenomena, involves the mass distribution of the 
wing and the damping properties of the structure). On the basis of 
past experience, certain qualitative conclusions can be drawn concerning 
this relation. 

As shown in the charts of this paper, the aileron reversal speed, 
or the speed at which a specified amount of control is retained, is 
lower, for highly sweptback wings than for unswept wings. Similarly, 
the divergence speed decreases rapidly as the angle of sweepforward 
increases. For a typical wing the values of q* at divergence and at 
reversal as obtained by the charts of this paper are shown as a func- 
tion of the sweep parameter k in figure 11. For unswept wings the 
dynamic pressure at flutter is usually within a certain range varying 
between a value lower than the dynamic pressure at reversal to a value 
higher than the dynamic pressure at divergence, depending on the 
geometric, structural, aerodynamic, and mass parameters of the given 
case, and varying even for a given case and a given speed range with 
altitude, because a change in air density may change the mode in which 
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the wing flutters. This wide range is indicated in figure 11 by- 
starting three flutter curves at k = 0; these curves do not necessarily 
constitute the upper and lower limits. 

If the variation of flutter speed with sweep angle is assumed for 
the purpose of illustration to he similar to that indicated in figure 17 
of reference 7, the three flutter curves of figure 11 are obtained. 

This figure must not be construed as presenting any quantitative infor- 
mation; to emphasize this point the wing is not identified. Even 
qualitatively the relation between the dynamic pressures at flutter, 
divergence, and reversal is subject to certain limitations because 
the flutter tests of referenced were performed 

(1) At subsonic speeds 

( 2 ) On models without ailerons 

- (3) On models without concentrated masses 

( 4 ) With models which fluttered in the classical two -degree -of - 
freedom mode 

There is reason to believe that sweptback wings with high aspect ratios 
flying at high altitudes may experience a possibly mild form of flutter 
in a single -degree -of- freedom mode, because a vertical motion necessarily 
implies vertical bending and, hence, in the case of a swept wing, a 
variation in the angle of attack. In general, the greater the number 
of degrees of freedom the more difficult it is to relate flutter to 
the static aeroelastic phenomena. 

However, at subsonic speeds and low or moderately high altitudes 
at least, the trend shown in figure II should be valid for wings with- 
out very large concentrated masses and with irreversible controls, 
which tend to minimize the possibility of aileron-coupled flutter. 
Consequently, if these wings are highly swept back they can be designed 
to meet lateral-control requirements with the likelihood that they will 
then be safe against flutter as well, provided conventional lateral- 
control devices are used. On the other hand, if these wings are 
unswept or even moderately swept back, they may have to be stiffened 
beyond the amount, required by static aeroelastic considerations. 

In any event the final design must be checked both for static 
aeroelastic effects and for flutter. In most cases the static aero- 
elastic effects can probably be calculated with sufficient accuracy by 
means of the charts and approximate formulas of reference 2 and the 
present paper. In some cases, however, particularly if these effects 
are in any way critical, a more refined method of analysis, such as 
that of references 1 and 8, may have to be used. 
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CONCLUDING REMARKS 


An approximate method based on charts and approximate formulas 
has been presented for estimating rapidly the aeroelastic effects on 
the lateral control of swept and unswept wings at subsonic and super- 
sonic speeds. The charts and approximate formulas presented in this 
paper together with those presented in NACA TN 2608 also serve to 
simplify design procedure in many instances because they can be used 
at the preliminary design stage to estimate the amount of additional 
material required to stiffen a wing which is strong enough and because 
th6y indicate that the best way of distributing this additional 
material is to locate most of it near the wing tip in most cases. 

For the purpose of making specific calculations, the limitations 
of the method of this paper are that they do not apply directly to 
wings with very low aspect ratio, with very large angles of sweep, with 
zero taper ratio, or with large sources of concentrated aerody nam ic 
forces . 

The charts and approximate formulas indicate that the control 
effectiveness of an airplane may be increased by varying some of the 
design parameters such as the ratio of torsional to bending stiffness 
and, if necessary, resorting to unconventional lateral-control devices. 
The charts also indicate that a wing which is strong enough is most 
likely to be affected by losses in lateral control due to wing flexi- 
bility if it is to operate at high dynamic pressures, if it is thin, if 
it has a large angle of sweepback, if it has an elastic-axis location 
relatively far forward on the chord or a location of the center of 
pressure due to aileron deflection far rearward on the chord, or if 
it is to operate at transonic or high-supersonic Mach numbers. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. , March 7> 1952 
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APPENDIX 


METHODS OF CALCULATIONS ON WHICH CHARTS ARE BASED 
The Aeroelastic Equations 


The assumptions made in the following analysis are the same as 
those made in reference 2: 

(1) Aerodynamic induction is taken into account by applying an 
over-all correction to strip theory 

(2) Aerodynamic and elastic forces are based upon the assumption 
of small deflections 

(3) The wing is clamped at the root perpendicular to a straight 
elastic axis, and all deformations are considered to be given by the 
elementary theories of bending and torsion about the elastic axis 

In addition it is assumed, as in reference 1 and elsewhere, that the 
angle between the aileron and the wing is constant along the span of 
the aileron. Then, for a wing with an outboard aileron, the force per 
unit width on sections perpendicular to the elastic axis is 

qcC]j a 

1 = + “S^/a) ( A1 > 


where _/• is a unit step function of s defined by 

Q, 





s t 
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where is the spanwise ordinate of the inboard end of the aileron. 

The running torque of this force about the elastic axis is 

') c2|! 1 C L 

£(a s - € 056 . 4 ) (A2) 


where e 


is the moment -arm ratio 



1 - 


The integration of these forces yields the accumulated torque and 
bending moment: 


T = 


q c r 2 eiC^ e 

144 





(A3) 



(A4) 


Combining these equations with the equations of elastic deformation 
presented in appendix A of reference 2 as 



results in two simultaneous equations of equilibrium: 



3 ^ 
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After equations (A5) and (A 6) have been solved, the rolling moment 
about the wing root may be found from the expression 

qs |c z 

— -T 7 — = T r sin A + M r cos A (A9) 


where the root twisting moment T p and the root bending moment Mj. 
are given by equations (A3) and (a 4) evaluated at s equal to zero. 
Then the rolling-moment ratio becomes 



1 + 


^(Bp + cBg) + B3 + eBl*. 


r eBc 


b 6 


(A10) 


where the functions B-^ 


to Bg are defined by 


B 1 ® Bg 



(Alla) 


b 3 + eB 4 



(Allb) 


b 5 



(Allc) 



(Alld) 
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and where the parameter d/k 


is defined by 


d 

k 


e l C r 


tan A 


(A12) 


Solution for Uniform Wings 

If the torsional stiffness, the bending stiffness, and the chord 
of the wing have constant values of (GJ) r , (El) r , and c r , respec- 
tively, along the wing span, the equations of equilibrium (A5) and (a 6 ) 
become 


<p * * cos A = 



cos A - T sin 


A) 



(A13) 


r»»» sin A 



cos A - T sin A) + 




(A14) 


where the differentiation denoted by the prime is with respect to 



Differentiating equation (M3) once with respect to £ and com- 
bining it with equation (klk) yields the single differential equation 
of equilibrium. 


a 3 '»» + q*a a * - qa s = a .56 q*e_7 ’ + qj (A15) 

a a 

subject to the following boundary conditions: 

o-s (!) = 0 (Al6a) 

a s *(0) = 0 (Al 6 b) 

a s * * (O) = -q* [ttg (O) - ?ag5_/ a ( 0 )J (Al6c) 
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The complete solution of equation (A15) can be readily obtained by 
means of Laplace transforms as 


a s(l) 

058 


= q*(l + 6) 


f 5 (n - f 5 d - ii) - £ ^ . 1} - ^ - 6i) f ^ (l) 


^(1) 


f 3d - ii) 

f^(i) 


f 3 (i) - f 3 (i - t ± ) - i a (D 


(A17) 


where the functions f 3 (|) and f^), as well as f^(|) which will 
be used subsequently, are defined in appendix A of reference 2 as 


f 3 (l) 

^(l) 


c l e “ 2 ^^ + cos 7l + — sin 7^ 

C^e -2 ^ + cos 7 1 + — sin 7|^ 


f^(4) = C^e~ 2 ^ + e^/cg cos 75 + — sin yl 

7 


where the constants of integration are defined in reference 2 in terms 
of the roots of the characteristic equation, -20 and 0 ± i7. These 
three functions are equal to zero when | < 0. 

The substitution of this solution into equations (All) yields the 
functions 


Bl(q * ,k) = ~rrr^ f ^ 1] - f ^ (1 - si) - ti + % (Ai8a) 


B 2 (q%k) = I (£3(1) - f 3 (l - hT] + q* 


1 f4(D 


k f*(l) 


[f 5 (l) - f 5 (l - ij] 


(Al8b) 
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B 3 (q*,k) 


f 3 (l ~ .iil f 5 (l) - f 5 (l - U) - Si + \ Si 2 + b 4 
f 3 (D 


(Al8c) 


B U (q*,k) = 4B (1) " f 3 (l “ + ^& (1) ' f ^ (l " 4i il 


_L . is 


f-s(l) 


k2 f 3 (l) 


(f 5 (i) - f 5 d - e ± )] - 1 ii (Ai8d) 


b 5 = -!i 


(Al8e) 


b 6 = 6i - \ ^ 


(Al8f) 


The value of q* at reversal is that value which makes c I§/ c I§q in 

equation (A10) equal to zero for. given values of the parameters k, e, 
and d. 

For a full-span aileron (li = 1), equations (Al8) become 


f 4 (l) 

Bi(q*,k) = 1 + Bg 

f 3 (!) 


(A19a) 


F 3 (1) - il + q* 

i mi) 

U> —J 

k f 3 d) 


f 5 (l) 


(A19b) 


B 3 (q*,k) 


f 5 (l) 

f 3 (l) 



(A19c) 
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B 4 (q*,k) = i(f 3 (l) - l] + i[? 4 (l) - 1 

1 ♦ 9 *[i - f 5 (1) > 5 (l) 
1 [fc2 f 3 (l)J 5 ' 

(A19d) 

b 5 = -1 


(A19e) 

B 6 = 2 


(Al9f) 


and for a "tip aileron" of very short span (l-^— >0), equations (18) 
become 

Bl(q*,k) = iq jf 3 (l) + q*[f 4 (l) - kf 5 (l^£it|il - l| + B 2 (A20a) 


B 2 (q*,k) = li 


q*f 5 (l) - q* f 4 (l) 

? Ml) 


(A20b) 


B 3 (q*,k) = li^f 4 (l) + q*[f 4 (l) - kf 5 (l)] 


^ 5 ( 1 ) 

MD 


If + B 4 (A20 c ) 


B 4 (q*,k) = sJ-iJfgU) - I] + q* 


1 £4^) 

k " ^ 


f 5 (l) 


(A20d) 


b 5 = -Si 


(A20e) 


B 6 = li 


(A20f) 
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Solution for Nonuniform Wings 

By means of strip theory applied at a finite number of points, 
equation (Al) may be written in matrix notation as 


{*} = -"j-jp l/ljfs} + “8 s { 4 }} 


(A21) 


and the expression for torque (A2) as 


« e l C I«e 


I41t 


[c 2 l 


{a 8 } - 


(A22) 


The matrix notation used in this analysis is the same as in reference 1. 
Equations (A3) and (A4), written in matrix notation, are 


W = Cl *e & ’] [(£) ] {w " i 1 ] ' 


(A23) 




(A24) 


where the integrating matrices &■] and [n •] are, respectively, 
for single and double integration from tip to root and are given in 
reference 1. The value of W at the matrix station nearest the 
discontinuity of _^ a can modified as in reference 1 in such a 
manner that premultiplication of by {i f -J and Ip'd yields 

the same area and moment about the wing root, respectively, as would 
be obtained by analytic integration. This modification may also be 
regarded as an attempt to round off the lift distribution near the 
inboard end of the aileron in a physically reasonable manner. 
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The combination of equations (A23) and (A24) with matrix expres- 
sions for and {V} - in terms of and •^m|’ (see appendix A 

of reference 2) yields the matrix- equilibrium equation, 

[111 - q*(A]j ja s ] = -q*a 5 5{p} (A25) 

which can be solved for The aeroelastic matrix [aJ is defined 

in reference 1 and the column matrix is defined by 

® (-6, 

After the rolling-moment ratio expressed by equation (A10) is set 
equal to zero, the condition for aileron reversal is 

f(Bi + eB 2 ) + B 3 + eBj^ + | eB^ + Bg = 0 

where the matrix equivalents of equations (All) become 


Bl+ te7a) 

(A2Tb > 

B 5 = -ll'lj (^)]{la} < A2 7c) 

B 6“ I ' J[r:| < A2 Td) 
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Therefore, the condition for aileron reversal expressed in matrix 
notation is 


£N| + L 11 *lj 1^:1 M = 065 | - L n ’lJ[^l {?a} 

(A28) 

Solving equation (A28) for a§5 and multiplying the resulting equation 
by {3} yields 


a 6 6 {p} = 



IN 

(4 

n 

+ 

[n*v 

l|4 

{4 

I' 1 

N 

0 

■ L 

. n, ij 

L4 

w 


(A29) 


The substitution of this expression for ag5 £ 3 } into the equilibrium 
equation (A25) yields 


{%} “ q *R M fs] 

where the aileron-reversal matrix M is defined by 


(A30) 


N - [Aj - 






!N(£) + H 1 *il L41 


(A31) 


The value of q* at aileron reversal can be found by the iteration of 
this aileron-reversal matrix. 
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Representation of Results by Approximate Formulas 

Approximate formulas, similar to those in reference 2, have been 
used to combine the results of the many computations indicated in this 
analysis; as in the case of the formulas presented in reference 2, 
those presented in this section are based on considerations of a semi- 
rigid wing. 

The functions B-^ to in equation (A10) have been found to be 
given quite accurately by the following approximate expressions: 


Bi = -Bg 




(A32a) 


b 2 = -b 6 


q*/q* D k 3 + K6 
1 _ 9lL 1 - K 2 k 
«*D 


(A32b) 



-B 6 


Q*/ (K 2 + K 7 )k 

i _ i - Kpk 

** D 


(A32c) 


b 4 = " b 6 


q*/ q* D k u 

i _ S*_ i - K 2 k 

q*D 


(A32d) 


. , , b 5 = ~ b S K 3 (A32e) 

where the factors K 4 to Ky depend on the aileron span, the taper ratio 
and the spanwise variation of the bending and torsion stiffnesses: The 
factor K 2 is independent of the aileron span, and the factor Ko is 
independent of the stiffness variation. ' 5 
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The following approximate formula may he obtained by the substitu- 
tion of equations (A32) into equation (A10): 



c *e 0 1 - i. 




(A33) 


where the value of q at aileron reversal is 


tod - K 2 k )(! - K 3 4 d) 

= 

1 + (K]| + KgK^d)? + K^d + Kg — d + Kyk 


(A34) 


With the use of the approximate formula presented in reference 2 , 


K, 


n* - 

q D “ 


1 - 


Kgk 


(A3 5) 


equation (A34) becomes 


K l( X ‘ K 3 f d ) 

q* R = — 

1 + . 0 % + KgKgdJe + K 5 d + Kg |-d + Kyk 


(A36) 


The accuracy of equation (A36) compared with the results calculated 
directly by the method of the preceding section is illustrated in 
figures 6 and 12. 

The results of reference 2 indicate that the damping- in-roll 
derivative may be expressed approximately by 




q D 


(A37) 
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and since the wing-tip helix angle due to roll is 


Pb c ?s 
. 2V = "% 

an approximate formula for rolling maneuverability is 


(A3 8) 


pb 

2V 



(A39) 


Figure 13 shows the approximate formulas (A3 3 ) and (A3 9 ) 
agreement with more accurately computed values. 


to be in good 


An approximate expression for the structural twist due to aileron 
deflection similar to the expressions for the structural twists due to 
geometrical angles of attack. given in reference 2 has been deduced from 
the results of the analysis in the preceding section: 


a s _ JiAd ( K 4 e + K 2 k ) f a " K 2 k 

a 5 5 1 - - 3 _ 1 - Kok 

% 


(a4o) 


where f & and Af r 


are functions of the spanwise coordinate s*, the 

The accuracy 


wing chord and stiffness variations, and the aileron span 
of equation (A4o) is indicated in figure 14. 
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Figure 1.- Definitions of geometric parameters. 
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Figure 3.- Concluded. 
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(b) Wings with moment-arm ratio 




59 



Figure 5* - Concluded. 
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(a) Stiffnesses given by constant -stress criterion for 
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Figure 7-- The angle -of-attack distribution functions f a and Af 
for aileron deflections, s*- = 0 . 8 . 
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Dimensionless distance along span, s* 


(l>) Stiffnesses related to those given by constant-stress criterion for 

wings of taper ratio 0 . 5 . 


Figure 7*- Continued. 
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Dimensionless distance along span, s 

(c) Stiffnesses proportional to cV 
Figure 7.- Concluded. 
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Figure 9.- Effect of aeroelastic action on some lateral-control properties 

of the example wing. 
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Figure 14. - Comparison of angle-of-attack ratios calculated by the matrix 
method of the appendix with those calculated by equation (a 40) for 
uniform wings, s*.j_ = 0.5; e =1.0. 
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